Single-nucleus and spatial transcriptomic profiling of human temporal cortex and white matter reveals novel associations with AD pathology.

Alzheimer's disease (AD) is a neurodegenerative disorder with complex pathological manifestations and is the leading cause of cognitive decline and dementia in elderly individuals. A major goal in AD research is to identify new therapeutic pathways by studying the molecular and cellular changes in the disease, either downstream or upstream of the pathological hallmarks. In this study, we present a comprehensive investigation of cellular heterogeneity from the temporal cortex region of 40 individuals, comprising healthy donors and individuals with differing tau and amyloid burden. Using single-nucleus transcriptome analysis of 430,271 nuclei from both gray and white matter of these individuals, we identified cell type-specific subclusters in both neuronal and glial cell types with varying degrees of association with AD pathology. In particular, these associations are present in layer specific glutamatergic (excitatory) neuronal types, along with GABAergic (inhibitory) neurons and glial subtypes. These associations were observed in early as well as late pathological progression. We extended this analysis by performing multiplexed in situ hybridization using the CARTANA platform, capturing 155 genes in 13 individuals with varying levels of tau pathology. By modeling the spatial distribution of these genes and their associations with the pathology, we not only replicated key findings from our snRNA data analysis, but also identified a set of cell type-specific genes that show selective enrichment or depletion near pathological inclusions. Together, our findings allow us to prioritize specific cell types and pathways for targeted interventions at various stages of pathological progression in AD.

Disease-specific selective vulnerability and neuroimmune pathways in dementia revealed by single cell genomics.

The development of successful therapeutics for dementias requires an understanding of their shared and distinct molecular features in the human brain. We performed single-nuclear RNAseq and ATACseq in Alzheimer disease (AD), Frontotemporal degeneration (FTD), and Progressive Supranuclear Palsy (PSP), analyzing 40 participants, yielding over 1.4M cells from three brain regions ranging in vulnerability and pathological burden. We identify 35 shared disease-associated cell types and 14 that are disease-specific, replicating those previously identified in AD. Disease - specific cell states represent molecular features of disease-specific glial-immune mechanisms and neuronal vulnerability in each disorder, layer 4/5 intra-telencephalic neurons in AD, layer 2/3 intra-telencephalic neurons in FTD, and layer 5/6 near-projection neurons in PSP. We infer intrinsic disease-associated gene regulatory networks, which we empirically validate by chromatin footprinting. We find that causal genetic risk acts in specific neuronal and glial cells that differ across disorders, primarily non-neuronal cells in AD and specific neuronal subtypes in FTD and PSP. These data illustrate the heterogeneous spectrum of glial and neuronal composition and gene expression alterations in different dementias and identify new therapeutic targets by revealing shared and disease-specific cell states.

A structured evaluation of cryopreservation in generating single-cell transcriptomes from cerebrospinal fluid.

Single-cell transcriptomics allows characterization of cerebrospinal fluid (CSF) cells at an unprecedented level. Here, we report a robust cryopreservation protocol adapted for the characterization of fragile CSF cells by single-cell RNA sequencing (RNA-seq) in moderate- to large-scale studies. Fresh CSF was collected from twenty-one participants at two independent sites. Each CSF sample was split into two fractions: one was processed fresh, while the second was cryopreserved for months and profiled after thawing. B and T cell receptor sequencing was also performed. Our comparison of fresh and cryopreserved data from the same individuals demonstrates highly efficient recovery of all known CSF cell types. We find no significant difference in cell type proportions and cellular transcriptomes between fresh and cryopreserved cells. Results were comparable at both sites and with different single-cell sequencing chemistries. Cryopreservation did not affect recovery of T and B cell clonotype diversity. Our CSF cell cryopreservation protocol provides an important alternative to fresh processing of fragile CSF cells.

Schizophrenia risk conferred by rare protein-truncating variants is conserved across diverse human populations.

Schizophrenia (SCZ) is a chronic mental illness and among the most debilitating conditions encountered in medical practice. A recent landmark SCZ study of the protein-coding regions of the genome identified a causal role for ten genes and a concentration of rare variant signals in evolutionarily constrained genes1. This recent study-and most other large-scale human genetics studies-was mainly composed of individuals of European (EUR) ancestry, and the generalizability of the findings in non-EUR populations remains unclear. To address this gap, we designed a custom sequencing panel of 161 genes selected based on the current knowledge of SCZ genetics and sequenced a new cohort of 11,580 SCZ cases and 10,555 controls of diverse ancestries. Replicating earlier work, we found that cases carried a significantly higher burden of rare protein-truncating variants (PTVs) among evolutionarily constrained genes (odds ratio = 1.48; P = 5.4 × 10-6). In meta-analyses with existing datasets totaling up to 35,828 cases and 107,877 controls, this excess burden was largely consistent across five ancestral populations. Two genes (SRRM2 and AKAP11) were newly implicated as SCZ risk genes, and one gene (PCLO) was identified as shared by individuals with SCZ and those with autism. Overall, our results lend robust support to the rare allelic spectrum of the genetic architecture of SCZ being conserved across diverse human populations.

Cell-type-specific cis-eQTLs in eight human brain cell types identify novel risk genes for psychiatric and neurological disorders.

To date, most expression quantitative trait loci (eQTL) studies, which investigate how genetic variants contribute to gene expression, have been performed in heterogeneous brain tissues rather than specific cell types. In this study, we performed an eQTL analysis using single-nuclei RNA sequencing from 192 individuals in eight brain cell types derived from the prefrontal cortex, temporal cortex and white matter. We identified 7,607 eGenes, a substantial fraction (46%, 3,537/7,607) of which show cell-type-specific effects, with strongest effects in microglia. Cell-type-level eQTLs affected more constrained genes and had larger effect sizes than tissue-level eQTLs. Integration of brain cell type eQTLs with genome-wide association studies (GWAS) revealed novel relationships between expression and disease risk for neuropsychiatric and neurodegenerative diseases. For most GWAS loci, a single gene co-localized in a single cell type, providing new clues into disease etiology. Our findings demonstrate substantial contrast in genetic regulation of gene expression among brain cell types and reveal potential mechanisms by which disease risk genes influence brain disorders.

Epigenomic priming of immune genes implicates oligodendroglia in multiple sclerosis susceptibility.

Multiple sclerosis (MS) is characterized by a targeted attack on oligodendroglia (OLG) and myelin by immune cells, which are thought to be the main drivers of MS susceptibility. We found that immune genes exhibit a primed chromatin state in single mouse and human OLG in a non-disease context, compatible with transitions to immune-competent states in MS. We identified BACH1 and STAT1 as transcription factors involved in immune gene regulation in oligodendrocyte precursor cells (OPCs). A subset of immune genes presents bivalency of H3K4me3/H3K27me3 in OPCs, with Polycomb inhibition leading to their increased activation upon interferon gamma (IFN-γ) treatment. Some MS susceptibility single-nucleotide polymorphisms (SNPs) overlap with these regulatory regions in mouse and human OLG. Treatment of mouse OPCs with IFN-γ leads to chromatin architecture remodeling at these loci and altered expression of interacting genes. Thus, the susceptibility for MS may involve OLG, which therefore constitutes novel targets for immunological-based therapies for MS.

CellMixS: quantifying and visualizing batch effects in single-cell RNA-seq data.

A key challenge in single-cell RNA-sequencing (scRNA-seq) data analysis is batch effects that can obscure the biological signal of interest. Although there are various tools and methods to correct for batch effects, their performance can vary. Therefore, it is important to understand how batch effects manifest to adjust for them. Here, we systematically explore batch effects across various scRNA-seq datasets according to magnitude, cell type specificity, and complexity. We developed a cell-specific mixing score (cms) that quantifies mixing of cells from multiple batches. By considering distance distributions, the score is able to detect local batch bias as well as differentiate between unbalanced batches and systematic differences between cells of the same cell type. We compare metrics in scRNA-seq data using real and synthetic datasets and whereas these metrics target the same question and are used interchangeably, we find differences in scalability, sensitivity, and ability to handle differentially abundant cell types. We find that cell-specific metrics outperform cell type-specific and global metrics and recommend them for both method benchmarks and batch exploration.

muscat detects subpopulation-specific state transitions from multi-sample multi-condition single-cell transcriptomics data.

Single-cell RNA sequencing (scRNA-seq) has become an empowering technology to profile the transcriptomes of individual cells on a large scale. Early analyses of differential expression have aimed at identifying differences between subpopulations to identify subpopulation markers. More generally, such methods compare expression levels across sets of cells, thus leading to cross-condition analyses. Given the emergence of replicated multi-condition scRNA-seq datasets, an area of increasing focus is making sample-level inferences, termed here as differential state analysis; however, it is not clear which statistical framework best handles this situation. Here, we surveyed methods to perform cross-condition differential state analyses, including cell-level mixed models and methods based on aggregated pseudobulk data. To evaluate method performance, we developed a flexible simulation that mimics multi-sample scRNA-seq data. We analyzed scRNA-seq data from mouse cortex cells to uncover subpopulation-specific responses to lipopolysaccharide treatment, and provide robust tools for multi-condition analysis within the muscat R package.

Tau Pathology Drives Dementia Risk-Associated Gene Networks toward Chronic Inflammatory States and Immunosuppression.

To understand how neural-immune-associated genes and pathways contribute to neurodegenerative disease pathophysiology, we performed a systematic functional genomic analysis in purified microglia and bulk tissue from mouse and human AD, FTD, and PSP. We uncover a complex temporal trajectory of microglial-immune pathways involving the type 1 interferon response associated with tau pathology in the early stages, followed by later signatures of partial immune suppression and, subsequently, the type 2 interferon response. We find that genetic risk for dementias shows disease-specific patterns of pathway enrichment. We identify drivers of two gene co-expression modules conserved from mouse to human, representing competing arms of microglial-immune activation (NAct) and suppression (NSupp) in neurodegeneration. We validate our findings by using chemogenetics, experimental perturbation data, and single-cell sequencing in post-mortem brains. Our results refine the understanding of stage- and disease-specific microglial responses, implicate microglial viral defense pathways in dementia pathophysiology, and highlight therapeutic windows.

Targeted Treatment of Individuals With Psychosis Carrying a Copy Number Variant Containing a Genomic Triplication of the Glycine Decarboxylase Gene.

BACKGROUND: The increased mutational burden for rare structural genomic variants in schizophrenia and other neurodevelopmental disorders has so far not yielded therapies targeting the biological effects of specific mutations. We identified two carriers (mother and son) of a triplication of the gene encoding glycine decarboxylase, GLDC, presumably resulting in reduced availability of the N-methyl-D-aspartate receptor coagonists glycine and D-serine and N-methyl-D-aspartate receptor hypofunction. Both carriers had a diagnosis of a psychotic disorder. METHODS: We carried out two double-blind, placebo-controlled clinical trials of N-methyl-D-aspartate receptor augmentation of psychotropic drug treatment in these two individuals. Glycine was used in the first clinical trial, and D-cycloserine was used in the second one. RESULTS: Glycine or D-cycloserine augmentation of psychotropic drug treatment each improved psychotic and mood symptoms in placebo-controlled trials. CONCLUSIONS: These results provide two independent proof-of-principle demonstrations of symptom relief by targeting a specific genotype and explicitly link an individual mutation to the pathophysiology of psychosis and treatment response.

Altered human oligodendrocyte heterogeneity in multiple sclerosis.

Oligodendrocyte pathology is increasingly implicated in neurodegenerative diseases as oligodendrocytes both myelinate and provide metabolic support to axons. In multiple sclerosis (MS), demyelination in the central nervous system thus leads to neurodegeneration, but the severity of MS between patients is very variable. Disability does not correlate well with the extent of demyelination1, which suggests that other factors contribute to this variability. One such factor may be oligodendrocyte heterogeneity. Not all oligodendrocytes are the same-those from the mouse spinal cord inherently produce longer myelin sheaths than those from the cortex2, and single-cell analysis of the mouse central nervous system identified further differences3,4. However, the extent of human oligodendrocyte heterogeneity and its possible contribution to MS pathology remain unknown. Here we performed single-nucleus RNA sequencing from white matter areas of post-mortem human brain from patients with MS and from unaffected controls. We identified subclusters of oligodendroglia in control human white matter, some with similarities to mouse, and defined new markers for these cell states. Notably, some subclusters were underrepresented in MS tissue, whereas others were more prevalent. These differences in mature oligodendrocyte subclusters may indicate different functional states of oligodendrocytes in MS lesions. We found similar changes in normal-appearing white matter, showing that MS is a more diffuse disease than its focal demyelination suggests. Our findings of an altered oligodendroglial heterogeneity in MS may be important for understanding disease progression and developing therapeutic approaches.

A junction coverage compatibility score to quantify the reliability of transcript abundance estimates and annotation catalogs.

Most methods for statistical analysis of RNA-seq data take a matrix of abundance estimates for some type of genomic features as their input, and consequently the quality of any obtained results is directly dependent on the quality of these abundances. Here, we present the junction coverage compatibility score, which provides a way to evaluate the reliability of transcript-level abundance estimates and the accuracy of transcript annotation catalogs. It works by comparing the observed number of reads spanning each annotated splice junction in a genomic region to the predicted number of junction-spanning reads, inferred from the estimated transcript abundances and the genomic coordinates of the corresponding annotated transcripts. We show that although most genes show good agreement between the observed and predicted junction coverages, there is a small set of genes that do not. Genes with poor agreement are found regardless of the method used to estimate transcript abundances, and the corresponding transcript abundances should be treated with care in any downstream analyses.

Marker chromosome genomic structure and temporal origin implicate a chromoanasynthesis event in a family with pleiotropic psychiatric phenotypes.

Small supernumerary marker chromosomes (sSMC) are chromosomal fragments difficult to characterize genomically. Here, we detail a proband with schizoaffective disorder and a mother with bipolar disorder with psychotic features who present with a marker chromosome that segregates with disease. We explored the architecture of this marker and investigated its temporal origin. Array comparative genomic hybridization (aCGH) analysis revealed three duplications and three triplications that spanned the short arm of chromosome 9, suggestive of a chromoanasynthesis-like event. Segregation of marker genotypes, phased using sSMC mosaicism in the mother, provided evidence that it was generated during a germline-level event in the proband's maternal grandmother. Whole-genome sequencing (WGS) was performed to resolve the structure and junctions of the chromosomal fragments, revealing further complexities. While structural variations have been previously associated with neuropsychiatric disorders and marker chromosomes, here we detail the precise architecture, human life-cycle genesis, and propose a DNA replicative/repair mechanism underlying formation.

Progress in Understanding and Treating SCN2A-Mediated Disorders.

Advances in gene discovery for neurodevelopmental disorders have identified SCN2A dysfunction as a leading cause of infantile seizures, autism spectrum disorder, and intellectual disability. SCN2A encodes the neuronal sodium channel NaV1.2. Functional assays demonstrate strong correlation between genotype and phenotype. This insight can help guide therapeutic decisions and raises the possibility that ligands that selectively enhance or diminish channel function may improve symptoms. The well-defined function of sodium channels makes SCN2A an important test case for investigating the neurobiology of neurodevelopmental disorders more generally. Here, we discuss the progress made, through the concerted efforts of a diverse group of academic and industry scientists as well as policy advocates, in understanding and treating SCN2A-related disorders.

Divergent Levels of Marker Chromosomes in an hiPSC-Based Model of Psychosis.

In the process of generating presumably clonal human induced pluripotent stem cells (hiPSCs) from two carriers of a complex structural rearrangement, each having a psychotic disorder, we also serendipitously generated isogenic non-carrier control hiPSCs, finding that the rearrangement occurs as an extrachromosomal marker (mar) element. All confirmed carrier hiPSCs and differentiated neural progenitor cell lines were found to be mosaic. We caution that mar elements may be difficult to functionally evaluate in hiPSC cultures using currently available methods, as it is difficult to distinguish cells with and without mar elements in live mosaic cultures.

Contribution of copy number variants to schizophrenia from a genome-wide study of 41,321 subjects.

Copy number variants (CNVs) have been strongly implicated in the genetic etiology of schizophrenia (SCZ). However, genome-wide investigation of the contribution of CNV to risk has been hampered by limited sample sizes. We sought to address this obstacle by applying a centralized analysis pipeline to a SCZ cohort of 21,094 cases and 20,227 controls. A global enrichment of CNV burden was observed in cases (odds ratio (OR) = 1.11, P = 5.7 × 10-15), which persisted after excluding loci implicated in previous studies (OR = 1.07, P = 1.7 × 10-6). CNV burden was enriched for genes associated with synaptic function (OR = 1.68, P = 2.8 × 10-11) and neurobehavioral phenotypes in mouse (OR = 1.18, P = 7.3 × 10-5). Genome-wide significant evidence was obtained for eight loci, including 1q21.1, 2p16.3 (NRXN1), 3q29, 7q11.2, 15q13.3, distal 16p11.2, proximal 16p11.2 and 22q11.2. Suggestive support was found for eight additional candidate susceptibility and protective loci, which consisted predominantly of CNVs mediated by nonallelic homologous recombination.

Translating genome-wide association findings into new therapeutics for psychiatry.

Genome-wide association studies (GWAS) in psychiatry, once they reach sufficient sample size and power, have been enormously successful. The Psychiatric Genomics Consortium (PGC) aims for mega-analyses with sample sizes that will grow to >1 million individuals in the next 5 years. This should lead to hundreds of new findings for common genetic variants across nine psychiatric disorders studied by the PGC. The new targets discovered by GWAS have the potential to restart largely stalled psychiatric drug development pipelines, and the translation of GWAS findings into the clinic is a key aim of the recently funded phase 3 of the PGC. This is not without considerable technical challenges. These approaches complement the other main aim of GWAS studies, risk prediction approaches for improving detection, differential diagnosis, and clinical trial design. This paper outlines the motivations, technical and analytical issues, and the plans for translating PGC phase 3 findings into new therapeutics.

Frequency and Complexity of De Novo Structural Mutation in Autism.

Genetic studies of autism spectrum disorder (ASD) have established that de novo duplications and deletions contribute to risk. However, ascertainment of structural variants (SVs) has been restricted by the coarse resolution of current approaches. By applying a custom pipeline for SV discovery, genotyping, and de novo assembly to genome sequencing of 235 subjects (71 affected individuals, 26 healthy siblings, and their parents), we compiled an atlas of 29,719 SV loci (5,213/genome), comprising 11 different classes. We found a high diversity of de novo mutations, the majority of which were undetectable by previous methods. In addition, we observed complex mutation clusters where combinations of de novo SVs, nucleotide substitutions, and indels occurred as a single event. We estimate a high rate of structural mutation in humans (20%) and propose that genetic risk for ASD is attributable to an elevated frequency of gene-disrupting de novo SVs, but not an elevated rate of genome rearrangement.

The Influence of Microdeletions and Microduplications of 16p11.2 on Global Transcription Profiles.

Copy number variants (CNVs) of a 600 kb region on 16p11.2 are associated with neurodevelopmental disorders and changes in brain volume. The authors hypothesize that abnormal brain development associated with this CNV can be attributed to changes in transcriptional regulation. The authors determined the effects of 16p11.2 dosage on gene expression by transcription profiling of lymphoblast cell lines derived from 6 microdeletion carriers, 15 microduplication carriers and 15 controls. Gene dosage had a significant influence on the transcript abundance of a majority (20/34) of genes within the CNV region. In addition, a limited number of genes were dysregulated in trans. Genes most strongly correlated with patient head circumference included SULT1A, KCTD13, and TMEM242. Given the modest effect of 16p11.2 copy number on global transcriptional regulation in lymphocytes, larger studies utilizing neuronal cell types may be needed in order to elucidate the signaling pathways that influence brain development in this genetic disorder.